Calculation of the expansion rate of the three-volume measure in 

high-energy heavy-ion collisions 

A. Dumitru 

Physics Department, Yale University 
P.O. Box 208124, New Haven, CT 06520, USA 

In ultrarelativistic heavy-ion collisions the local three-volume measure is 
expanding in the longitudinal and transverse directions. This is similar to 
the Hubble-expansion of the universe, except that the former is not locally 
isotropic. As an example the expansion rate is calculated assuming that the 
energy-momentum tensor in the central region is that of an ideal fluid, un- 
dergoing Bjorken flow in longitudinal direction, and with initial conditions as 
expected for BNL-RHIC energy. While the longitudinal expansion of three- 
volume is independent of the energy density of the fluid, in case of 3+1 di- 
mensional expansion the form of the hydrodynamical solution (rarefaction 
CN ' wave or deflagration shock) affects the three-volume expansion rate on the 

■ hadronization hypersurface. As a consequence the average expansion rate on 
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that surface depends on the transverse size of the system. This may reflect in 
an impact-parameter dependence of the formation probability of light nuclei 



in , 

and of the freeze-out temperature of the strong interactions in the system. 

Q\ 

In high-energy hadron-hadron or nucleus-nucleus collisions secondary particles are sup- 
posedly produced on a so-called proper-time hyperbola r = \/t 2 — z 2 =const., with their 
number and energy density distributions independent of space-time rapidity r\ = | log |^ [T|]. 
(z and t denote the longitudinal coordinate and the time measured in the global rest-frame 
of the reaction, respectively.) This should be a reasonable (qualitative) approximation at 
^ ■ high energies and around midrapidity. In other words, secondary particle production oc- 
curs such that invariance under longitudinal Lorentz boosts, i.e. rapidity shifts, is obeyed. 
E.g., this holds true at asymptotically high energies for classical non-Abelian Yang-Mills 
Bremsstrahlung emitted by sources of color-charge on the light-cone @, if recoil is ne- 
glected. The effect is also built in effective string-models for particle production at high 
energies: a hadron emerging from the string-decay with rapidity y p gets on mass-shell in the 
global rest-frame at time t = TfCoshy p , where Ty is the formation time of the hadron in its 
rest-frame 0. For a discussion within the parton model see ref. [|J. 

The subsequent dynamical evolution preserves the invariance of the bulk properties, e.g. 
the energy density distribution, under longitudinal Lorentz boosts. It is therefore convenient 
to switch from (t, x, y, z) to new coordinates, via |l] 

£ = tcos1i?7 , z = rsinh?7 , x = r±cos<j) , y = r±sin<j) . (1) 

The Minkowski line element in terms of the new coordinates is ds 2 = dr 2 — r 2 di] 2 — dr\ — 
r^d^ 2 , i.e. g^ v = diag(l, — r 2 , — 1, — r\). In the following, vectors will be written in the 
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basis were the components 0, 1, 2, 3 correspond to the r, rj, r±, <p direction, respectively. 
Furthermore, rotational symmetry around the longitudinal direction is assumed. 

Given some boundary conditions on the hypersurface r = r , the evolution within the 
forward light-cone is governed by Q 

d ■ T = , d ■ N B = . (2) 

In the ideal fluid approximation (natural units are employed, h = c = k B = 1) 

T» u = (e + p)u»u v - V cT , N£ = n B u» . (3) 

The four-velocity of the locally comoving frame is normalized to u ■ u — 1, and is given by 
m m = 7_l(1, 0, v±, 0), where 7J 2 = 1 — v \ 0-0. The projections of the equations for the 
energy-momentum tensor parallel and orthogonal to u, and the continuity equation for the 
net baryon current yield 

u ■ de = — (e + p) d ■ u , (4) 
d vP = , (5) 
u-dn B = —n B d-u , (6) 

with 

u ■ d = 7 ± (d T + v ± d ± ) . (7) 
The expansion scalar d ■ u is given by 



g \ ' Qt T Tj_ Or 



Thus, for purely longitudinal expansion (u± — 0, 7j_ = 1), we simply have d ■ u = 1/r. This 
is obvious because in this case V = ttR^t is the three-volume on a r = const, hypersurface 
corresponding to a length in space-time rapidity of Arj = 1. In other words, the Hubble 
constant in longitudinal direction is 

H(r) = - , (9) 

T 

for r > r and rj_ < R±. H(t ~ 0.5 — 1 fm) is roughly 10 40 times bigger than the present rate 
of expansion of the universe, and about 10 18 times the Hubble constant at the cosmological 
hadronization phase transition. Note that for purely longitudinal expansion d ■ u (or H) is 
independent of e and n B , i.e. the fluid does not influence the expansion rate. This is due to 
eq. @: free (non-accelerated) flow always means 77 = const., independently of the equation 
of state p(e, n B ). 

As a side-remark, note that for spherically symmetric, three-dimensional boost-invariant 
expansion, v = f/t. One finds d ■ u = 3/r, where ([[]) is now replaced by t — r cosh rj, 
r = rsmhi], that is r 2 = t 2 — r 2 and rj = |log^. Again, the expansion scalar d ■ u is 
independent of the equation of state. 
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One also observes from eqs. (f|,|6|) that the continuity equations (0) can not be extrapolated 
to r = because the three- volume vanishes (the Jacobian \ (d(t, z)/d(r, r])\ = r of ([!]) is zero). 
Nevertheless, the transformation from (t, z) to (r, r/) is not useless, because the classical 
description breaks down for r — > anyway. In such reactions the {p±) of produced quarks 
and gluons are on the order of 1 GeV, and thus the uncertainty relation sets a time-scale of 
~ 1/ (p±) = 0.2 fm, below which this approach must fail. 

If transverse expansion is superimposed on longitudinally boost-invariant expansion, d ■ u 
of course depends on the transverse flow velocity, cf. eq. (||), which by itself depends on 
transverse pressure gradients (caused by energy density and/or baryon density gradients). 
Thus, the evolution of e, n B , u, and d ■ u is coupled via eqs. (§,§J). Consider, in particular, a 
hypersurface = (r, i], r_i_, (p) in space-time (e.g. a surface of constant time, or temperature, 
or the hadronization hypersurface A = 0, where A denotes the local fraction of quark-gluon 
phase). In parametric representation, is a function of three parameters ||. In our case, 
due to the symmetry under rotations around and Lorentz-boosts along the beam axis, two 
of these parameters can simply be identified with rj and <ft, while r and r± depend only on 
the third parameter, call it (. Thus, ( G [0, 1] parametrizes the hypersurface in the planes 
of fixed rj and (counter clock- wise). Then, the normal is[] 

du a du^ do" 1 ( dr± 8t \ 

da M = — — — dCdr/d0 = r ± r I -— , 0, — , 1 dCdr/d0 (10) 

The three-volume measure on the hypersurface is 

(Qt dr± \ 

E.g., on constant-r hypersurfaces and for purely longitudinal flow, d^ = Td 2 r±dr]. Note 
that d ■ u is simply the rate of expansion of the three- volume measure (|TT|), since 



u -ddV = dVd-u . (12) 

Eq. (|TJ) can be verified by an explicit calculation using the relations (|7|J^,p!l|) . However, the 
following proof is simpler and more general (it does not assume longitudinal boost-invariance 
and cylindrical symmetry). Note that the total net baryon number 

B = J da^ = J dVn B (13) 

is a constant, and thus u ■ dB = 0. Therefore, 

J (u ■ ddV) n B = u-d J dVn B - J dVu ■ dn B = J dVn B (d ■ u) . (14) 



x In our basis contains a factor ^J—g = r±r. 
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In the second step, the continuity equation (|) has been used. This must be true for any 
arbitrary function rag, and therefore the identity ( |12"D must hold. 

To discuss a specific example, equations (0) have been solved numerically employing 
the finite-difference scheme RHLLE, as described and tested in ||10| . The initial conditions 
were chosen as might be appropriate for collisions of heavy nuclei at BNL-RHIC energy, 
y/s = 200 GeV per incident nucleon pair. The main goal [[[]]] of these experiments is to 
repeat, in the laboratory, the QCD phase transition that occured at some stage in the 
universe. As already mentioned above, the difference is that the "Hubble-constant" is much 
larger in high-energy heavy-ion collisions. The energy density is expected HTT| , |T2ll to be 



significantly higher than in lower-energy reactions, such that the expansion effect, and the 
influence of the QCD hadronization phase transition, should be more prominent. 

In particular, on the r = r hypersurface, an entropy per net baryon of s/ns = 200 is 
assumed; see || for the resulting single-particle transverse momentum spectra as well as 
average transverse momenta and velocities of various hadrons. For simplicity, s and ns are 
assumed to be homogeneously distributed on the r = r hypersurface, independent of rj and 
0, and proportional to a step function, 6 (Rj_ — r±). Except close to the light-cone, where 
the flow velocities are largest but where the fluid is very dilute, the bulk dynamics is not 
very sensitive to the precise initial profile |J. 

Estimates for tq at RHIC energy span the range 0.2 — 1 fm [fT2" |. Here, the average value 



of To = 0.6 fm is employed. For the given initial conditions, and for the MIT bag-model 
equation of state (EoS) described below, the initial temperature is T(r ) ~ 300 MeV. 

To close the system of equations (fj) one has to specify an EoS, which determines the 
function p(e, n^). In the low-temperature region a gas of non-interacting hadrons that in- 
cludes all known |13[ strange?] and non-strange hadrons and hadronic resonances up to a 
mass of 2 GeV is assumed. At high temperatures, the EoS is that of an ideal gas of quarks, 
antiquarks (with masses m u = = 0, m s = 150 MeV), and gluons. In addition, a bag 



term [14] +B in the energy density, and —B in the pressure is included. This is the same 
as adding a term Bg^ v to the energy-momentum tensor of the quark-gluon fluid. Math- 
ematically, the bag term resembles the cosmological-constant term which is introduced in 
some cosmological models [||]. However, since it is added to the energy-momentum tensor of 
the quark-gluon fluid only, and not to that of the hadronic fluid, it changes the form of the 
hydrodynamical solution (deflagrations can occur, see below). 

Finally, the phase coexistance region is obtained from Gibbs' conditions of phase equi- 
librium. Thus, the EoS exhibits a first-order phase transition. The value of B is determined 



by Tc\ for Tc = 160 MeV, which is roughly in accord with recent lattice-QCD results [|15 
one finds B = 380 MeV/fm 3 . 



2 The strangeness chemical potential is determined by the requirement that the net strangeness 
density ns, and thus the strangeness current Ng = nsu^, vanish everywhere. 
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FIG. 1. Hypersurfaces corresponding to A = 1 (boundary between QGP and mixed phase), A = 
(boundary between mixed phase and pure hadron phase, i.e. the hadronization hypersurface) , and 
the T = 140 MeV isotherm; for t /R± = 1/10. 

Fig. |l| shows the hypersurfaces where the pure QGP and the mixed phase end, respec- 
tively, as well as the T = 140 MeV isotherm. One observes a rarefaction wave developping at 
rj_ = R±, r = tq. It accelerates a part of the quark-gluon fluid in transverse direction, before 
the boundary between quark-gluon plasma and mixed phase is reached (A = 1). Those parts 
of the system which have not been affected by transverse expansion enter the mixed phase 
at the same time, which is the horizontal part of the A = 1 hypersurface at r/R± ~ 0.6. 
The expansion rate within those space-time regions is simply 1/r, as discussed above. 

During the transition to hadronic matter, entropy is converted from the quark-gluon to 
the hadronic phase at a rate which of course depends on d ■ u. For example, assume for 
simplicity a net baryon free fluid, and that shock-solutions do not occur, such that entropy 
is conserved. Then, the continuity equation for the entropy four-flow, d ■ (su) = 0, yields 



u 



d\ = -(d-u)[\ + 



M 



(15) 



s H , s Q denote the entropy densities of the hadronic and quark-gluon fluids at T = T c . 
A is the local fraction of quarks and gluons within the mixed phase, such that the total 
entropy density is s = \s Q + (1 — \)s H . Eq. ( |l~5"D shows that the rate of adiabatic conversion 
of quark-gluon fluid into hadronic fluid is governed by the three-volume expansion rate 
d ■ u, which is larger than 1/r in the space-time regions where transverse expansion is 
active. Therefore, hadronization is faster close to the surface, r±_ ~ R±, than in the interior, 
cf. Fig. |l|. In particular, the local expansion rate of the hadronization volume (A = 
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hypersurface) determines whether the emerging hadrons can maintain local equilibrium or 
not. 

Fig. |2] shows the expansion factor £ = rd ■ u as a function of r and r±. d-uis multiplied 
by t to have £ = 1 in space-time regions where the expansion is purely longitudinal. In the 
front-left corner one can see the expansion induced by the rarefaction in the quark-gluon 
fluid. Once the mixed phase is reached, the rarefaction wave "stalls". 

This is due to the very small velocity of sound in the phase coexistance region. For 
recent discussions of the consequences of this effect in cosmology see e.g. [IB]. In heavy- ion 
collisions, it leads to "stall" of the transverse expansion 0-0. Fig. shows that therefore £ 
is nearly time independent until hadronization is completed and the purely hadronic phase 
is reached (A = 0): the curves of constant expansion rate (right panel of Fig. |2|) at « 0.6 
and 0.7 are nearly parallel to the time-axis. 



r/R T 





FIG. 2. Left: the expansion factor £ = rd-u as a function of r and r±. Right: curves of constant 
£, labeled by the value of £ to which they correspond to; the thick lines depict the hadronization 
hypersurface, A = 0, and the T = 140 MeV isotherm, respectively. Both for tq/R± = 1/10. 

The hadronization hypersurface (A = —0, i.e. infinitesimally smaller than 0) and the 
T = 140 MeV isotherm from Fig. [I] are shown again in Fig. |2|. Clearly, the expansion rate 
on these hypersurfaces is not constant, i.e., they are not hypersurfaces of homogeneity. This 
is in contrast to an isotropically expanding homogeneous universe ||17|| . Thus, in very high- 
energy heavy-ion collisions it does not seem natural to assume decoupling on a hypersurface 
of constant temperature, as already pointed out in Jl8| . 



For the equation of state with a first-order phase transition, hadronization is accom- 
plished via a deflagration shock [fF||3l, |10||T9l| . (For the geometry at hand, the shock front is 
actually curved, as can be seen in Fig. ||]). The final-state of the shock is the steady-state 
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corresponding to the Chapman- Jouget point [|T(],[2(J (maximum entropy production). The 
hadronization point A = 0, however, is in general located on the shock (and has higher 
pressure than at the CJ-point). On the shock front £ —>■ oo for an ideal fluid. In reality, a 
finite viscosity will smear out the shock front slightly; in the figure, £ remains finite because 
in the numerical solution the derivatives in eq. (||) are replaced by finite differences. Behind 
the shock front the expansion proceeds again via simple rarefaction waves, and £ is finite. 

The large expansion rate of the hadronization three-volume could open the interest- 
ing possibility that in the vicinity of that hypersurface the fluid breaks up into smaller 
droplets |2T[| which decouple from each other. This can lead to rapidity fluctuations |pljp2 | 
or even the formation of droplets of disoriented chiral condensate |23| . However, this scenario 
will not be discussed here in detail. 




U. I KJ.Z. U.J U.t U.J U.U KJ.I W.O U.S I 
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FIG. 3. Curves of constant £ (for tq/R± = 1/25, i.e. larger system than in Fig. |2|) labeled by 
the value of £ to which they correspond to; the thick lines depict the hadronization hypersurface, 
A = 0, and the T = 140 MeV isotherm, respectively. 

Clearly, those three-volume elements of the hadronization hypersurface that lie on the 
shock front will expand rather rapidly, even if the shock-width were slightly smeared out by 
a finite viscosity. It is interesting to note that for the smaller initial radius one obtains £ > 1 
almost over the entire hadronization hypersurface (and behind it), while for large R±, cf. 
Fig. [| £ = 1 over a large part of the hadronization hypersurface. For very large systems the 
expansion scalar d ■ u approaches the "trivial", EoS independent form d ■ u — 1/r. 

Of course, the question arises how the large expansion rate affects the observed spectra 
of hadrons. The most pronounced effect would be expected in space-time regions where the 
expansion scalar is large and the fluid is already rather dilute. E.g., the decoupling from local 
thermal equilibrium should be influenced [p!8|,p4],p5| : a large d ■ u hinders scattering between 
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the particles of the fluid. Therefore, the produced hadrons should decouple at higher T as 
compared to the case of small expansion rate, cf. also the discussion in f|i8p4| . In other 
words, the strong interaction freeze-out hypersurface should approach the hadronization 
hypersurface as beam energy increases (as observed in p6[) and as the initial radius decreases. 

If indeed the strong interactions cease at higher temperature as compared to a slowly 
expanding three- volume measure, the average phase- space density of frozen-out pions or 
kaons should also increase. This quantity can be estimated from measurements of the two- 
particle correlation function | 27|| . Moreover, the synthesis of light nuclei (and anti-nuclei) 
might be suppressed if the three-volume measure expands rapidly, because it becomes more 
difficult to coalesce (anti-)nucleons into clusters. Also, unlike in a static fireball in global 
thermodynamical equilibrium one may expect that the decoupling conditions of various 
hadron species differ: a large volume expansion rate should lead to earlier decoupling of 
hadron species with small cross section . Finally, if the QCD chiral phase transition is of 
first-order only at larger baryon-chemical potential \xb @, but a smooth cross-over at small 
fiB, the hydrodynamical solution (and thus the expansion rate) might look very different on 
the two sides of the critical point. In the region where the phase transition is first order, a 
shock can develop and lead to a rather large expansion rate. In the cross-over region, on the 
other hand, only rarefaction waves but no shocks will form, and the average expansion rate 
at chiral symmetry restoration should be smaller. As discussed above, this difference might 
reflect in the freeze-out properties of the strongly interacting system. 

To analyze these points in detail, one has to supplement the fluid-dynamical solution on 
the hadronization hypersurface with a more detailed kinetic treatment of the hadronic stage, 
which explicitly accounts for the various (elastic and inelastic) elementary hadron-hadron 



scattering processes [p6 |. One might then be able to study the effect of the three- volume 



expansion on the evolution and freeze-out properties of the hadron fluid. Those calculations 
found that some hadrons are indeed emitted directly from the hadronization hypersurface 
without scattering any further, and that the large expansion rate almost "freezes" the chem- 
ical composition of the hadron fluid. More detailed studies, e.g. with varying t /Rj_, will be 
performed in the future. 

Similar studies could be done in the laboratory by varying the initial transverse size. One 
should keep in mind, however, that the energy density in the central region decreases with 
mass number, which could partly counterbalance the effect. 

In summary, transverse expansion couples the scalar d ■ u, where u is the four-velocity 
of the locally comoving frame, and which can be interpreted as the expansion rate of the 
local three- volume measure, to the properties of the fluid (energy and baryon density and the 
EoS). If hadronization proceeds via a shock wave, the expansion rate can become particularly 
large on the hadronization hypersurface. The hadrons produced on such a shock may decou- 
ple immediately. The average expansion rate on the hadronization hypersurface depends on 
the size of the system (at fixed initial energy density). This may reflect in a system-size de- 
pendence of the coalescence probability of light (anti-)nuclei and of the freeze-out properties 
of the strongly interacting system (average phase space density and temperature). 
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